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The quantum formalism permits one to discriminate sometimes between any set of linearly- 
independent pure states with certainty. We obtain the maximum probability with which a set of 
equally-likely, symmetric, linearly-independent states can be discriminated. The form of this bound 
is examined for symmetric coherent states of a harmonic oscillator or field mode. 



OO . I. INTRODUCTION 

G\ ■ 

It is possible to manipulate the state of a quantum system in far more interesting ways than can be achieved by 
carrying out unitary operations and von Neumann measurements on the system of interest alone. Consideration of 
the effects of interactions with other systems has led to the development of the quantum operations formalism [[l], 
which allows any completely positive, trace-preserving map to represent, in principle, a realisable transformation of 
the density operator p. One particularly interesting type of operation is a probablistic operation. This is an operation 
which, with some probability less than 1, will transform the state of the system in a manner which cannot be brought 
about by any deterministic process. Although such operations generally have a non-zero failure probability, one 
^ generally knows whether or not the desired transformation has taken place. 

An important class of probablistic operations are those which allow one to discriminate unambiguously between 
non-orthogonal states, that is, with zero probability of error. When carried out on a quantum system prepared 
in one of the non-orthogonal states \ipj), such an operation will, with some probability, transform the state into a 
corresponding member of an orthonormal set \(f>j). The latter states can be discriminated without error using a simple 
OO ■ von Neumann measurement. Although such an operation cannot have unit probability of success, we can always tell 
0^ ] whether or not the desired transformation has taken place. When the attempt fails, we obtain an inconclusive result. 

The subject of unambiguous state discrimination was pioneered a decade ago by Ivanovic Q , Dieks J3| and Peres [Q , 
O |' and has recently undergone interesting further developments. While earlier work concentrated only on the problem of 
discriminating between two states, the problem of discriminating between multiple states has since been addressed. In 
particular, one of us Q has shown that the necessary and sufficient condition for a set of states \ipj) to be amenable to 
unambiguous state discrimination is that they must be linearly-independent. More recently, the problem of discrimi- 
^ ( - nating between three states has been examined in detail by Peres and Terno ||. We have shown that unambiguous 
discrimination is intimately related to other well-known types of probablistic operation, such as entanglement concen- 
tration [^U^] and exact cloning [|| . The latter connection has also been examined by Guan and Duo ^] . Unambiguous 
discrimination between two non orthogonal states has been demonstrated in the laboratory by Huttner et al . In 
this experiment, weak pulses of light were prepared in non-orthogonal polarisation states, a fraction of which were 
converted into orthogonal ones by a loss mechanism. 

It is clearly of interest to find the optimum strategy for discriminating unambiguously between a set of known 
states, that is, to determine the maximum probability of obtaining a definite result. As with problems in conventional 
quantum detection theory [jnj, where the aim is to find the absolute maximum of the discrimination probability or 
mutual information for a given source, few analytic solutions for optimum strategies and their respective figures of 
merit are known. The complete solution for unambiguous discrimination between two states with arbitrary a priori 
probabilities has been found by Jaeger and Shimony GJ. Peres and Terno || have explored the geometry of the 
optimisation problem for three states, and obtained useful insight into the general N state case. As yet, however, no 
analytical solutions have been found for more than two states. Such a solution is given in this paper. We determine 
the maximum probability with which TV symmetric states can be unambiguously discriminated, assuming they have 
equal prior probabilities. We then apply our result to examine the maximum probability of discriminating between 
N symmetric coherent states \aj). 
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II. UNAMBIGUOUS STATE DISCRIMINATION 



Consider a quantum system prepared in one of N pure quantum states \ipj), where j = 0, N — 1. These states 
span an iV-dimensional Hilbert space TL. If the states are non-orthogonal, no quantum operation can deterministically 
discriminate between them. It is, however, possible to devise a strategy which, with some probability, will reveal the 
state with zero error probability. To see how this may be done, it is convenient to employ the Kraus representation 
of quantum operations ijp. Each of the possible, distinguishable outcomes of an operation, which are labelled by the 
index /i, is associated with a linear transformation operator A^. These form a resolution of the identity 

£4^ = 1 (2.1) 

If the system is prepared with the initial density operator p, the probability P^ of the pth outcome is Tr pAj^A^. 
The final density operator corresponding to this result is A^pA^/ P M . The state discrimination operation will have 
N + 1 distinct outcomes, corresponding to detection of each of the states, and an additional answer which gives no 
information about the state. This is the inconclusive result. The operator which corresponds to the detection of the 
state \ipj) is Aj, where j = 0, N — 1, and we let Ap be the operator which leads to a failure of the discrimination 
attempt. Clearly, we have 

A F A F + J2^j = ^ ( 2 - 2 ) 



The zero-errors condition takes the form 



{^ j \A] l A r \i> i )=P s S ir , (2.3) 



where Pj is the conditional probability, given that the system was prepared in the state \ipj), that this state will be 
identified. This zero-errors condition can only be met if the \tpj) are linearly- independent, and we find that the Aj 
have the form ||): 



pl/2 



where the \<j>j) form an orthonormal basis for TL. Here we have also introduced the reciprocal states The 
reciprocal state is defined as that which lies in TL and is orthogonal to all \ipj>) for jy^j'. The set of states \tpj 



is simply the TV-dimensional complex generalisation of the set of (normalised) reciprocal vectors in crystallography 
|l3f with respect to the unit cell basis vectors, these being the \ipj). A complete set of reciprocal states exists if, and 
onlv^ if , the are linearly- independent. The reciprocal states are also necessarily linearly-independent, as is shown 

Given that the states have a priori probabilities rjj, the total probability of correctly identifying the state is 

3 3 

It convenient to proceed using the language of Positive Operator- Valued Measures (POVMs) Qj. The measurement 
can be expressed as an N + 1 element POVM operation by defining the positive Hermitian operators Edj — AjAj 
and Ep = A^ f Af- It is also useful to define Ed = J2j AjAj. The discrimination probability Pd is constrained by 

the fact that Ep must be positive. This, together with the decomposition of the identity, Ep + Ed = 1, means that 
none of the eigenvalues of Ed may exceed unity. It has been shown that the optimum measurement corresponds to 
the maximum eigenvalue of Ed being equal to 1 B . 



III. MAXIMUM DISCRIMINATION PROBABILITY FOR SYMMETRIC STATES 

In this section we derive the maximum attainable value of the unambiguous discrimination probability Pd for 
symmetric states with equal a priori probabilities. A set of quantum states spanning a Hilbert space TL is 
symmetric ]l4|] if there exists a unitary transformation U on TL such that 
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\^) = U\^. X ) = W\%), (3.1) 
|Vto) = %iv-i>, (3.2) 
U N = i. (3.3) 

Eq. (3.3) follows from the fact that any state in Tt can be written as a superposition of the and from U N \ipj) — 
\ipj). Such states are also said to be covariant with respect to U and have been found to have a preferential status 
with regard to problems in quantum detection theory. In conventional quantum detection theory, where the aim 
is to maximise the discrimination probability or mutual information for all possible measurements and not just the 
subset defined by the no-errors constraint, the maximum discrimination probability for a set of symmetric states with 
equal a priori probabilities can be obtained exactly jll 14|. The optimum strategy uses the so-called 'square root' 



fj"T| , [i"4| or 'pretty good' |l5| measurement. Considerable progress has also been made towards maximising the mutual 
information for these states |]l6| , in particular in connection with symmetric quantum channels fll^| . 

As has been shown in ||, unambiguous discrimination between linearly- independent symmetric states arises nat- 
urally in connection with entanglement concentration, that is, transforming a fraction of an ensemble of systems all 
prepared in the same imperfectly entangled state into a maximally entangled state, using only local operations and 
classical communication. The protocol given there, a generalisation of the 'Procrustean' technique due to Bennett et al 
p8[ , will maximise the entanglement of a pair of subsystems with probability equal to the probability of discriminating 
between a certain set of linearly-independent symmetric states. 

Prior to solving for the maximum unambiguous discrimination probability for equally-probable linearly-independent 
symmetric states, we shall obtain a representation of them which simplifies our analysis. We use the fact that the 
operator U can be expanded as follows: 

N-l 

U= 5>**bfe><7*|. (3.4) 

fc=0 

where ("fk\jk') = 5kk' ■ The real angles <pk may be taken to lie in half-open interval [0,2tt). It follows from Eq. (3.3) 
that 



(3.5) 



where fk is an integer satisfying 0</fe<n — 1. It is convenient to arrange the fk in increasing order, so that 

for k > kl . Clearly, we can expand \ipo) as J2k Ck \^ k )' f° r some c fc satisfying J2k l Cfc l 2 = 1- Together with Eqs. (3.1) 

and (3.4), this leads to 



N-l 

- <*Mk 



fc=0 

Note that the linear independence of the l^) implies that all of the c k are non-zero. Linear independence means that 
no superposition of the \tpj) can vanish, so consider 

JV-l JV-l 

£ e^%> = J2 CkSrhhk). (3.7) 

3=0 k=0 

If the \if)j) are linearly-independent, this must be non-zero for all r — 0, ...N — 1. Therefore, the fk must take every 
value in this range of integers. As we have arranged these integers in increasing order, we find that fk is simply equal 
to fc, so that linearly-independent symmetric states necessarily have the form: 

N-l 

|V»j)= ]Tc fe e^| 7fc >. (3.8) 

fe=0 

That having this form is also a sufficient condition for linear independence is proven in ||. For these states, the 
corresponding reciprocal states are given by 

N-l 

l#) = ;r 1/2 £ c *"V^l7r> (3.9) 

r=0 



3 



where Z = \ C A 2 ■ Note that the \ipj) are also symmetric, with respect to the same transformation U as the \ipj). 
For symmetric states, the operator Ed has the explicit form 

^ = ^ H P^c-^e^^ 1 ^)^. (3.10) 

j,r,r' 

Let us denote by E^f 4 an operator of this form which gives the maximum value of Pjj. We do not assume this 
operator to be unique, that is, we do not assume the optimum Pj to be unique. We can however, show that there 
exists an optimal operator E° D pt which possesses the symmetry 

= UE^UK (3.11) 

We prove this by contradiction, by first supposing that no £^f satisfies Eq. (3.11). Consider now any operator Ed 
of the form (3.10). If Ed corresponds to the maximum value of Pd, then its maximum eigenvalue, which we denote 
by X + (Ed), is equal to 1. We then define 

E$ = U l E D U n . (3.12) 

These operators clearly have the same eigenvalues as Ed and give the same value of Pd- In fact, P5 can be obtained 
from Ed by cycling the probabilities Pj . Writing explicitly the dependence of these operators on the Pj , we see that 
En{Pj) = Ed(P 3 -i), where P j±N = Pj. Consider now the operator 

i N — l 

fr-jE^-yEniWW < 3 - 13 > 

The second equation here is true if the a priori probabilities rjj are all equal to 1/N, which we take to be the case. The 
operator E D ve is invariant under the similarity transformation E D ve ^UE D ve W and gives the same transformation 
probability as Ep. Forming E D ve from Ed amounts to replacing all of the Pj by their average value Pd, the quantity 
we wish to maximise. However, its maximum eigenvalue \+(E D ve ) satisfies 

A + (E% e )<^J2 X +^) = l- (3-14) 
i 

This is a consequence of the fact the the maximum eigenvalue is convex on the space of Hermitian operators on Ti. 
A simple proof of this is given in the Appendix. Let us define E D = E D ve /\ + (E D ve ). The maximum eigenvalue 
of this operator is 1, so it is physically admissable. However, the success probability for this operator is P D = 
Pd / A+{E D ve )>PD, from (3.14). Therefore, we can obtain from any Ed another operator E D which is invariant under 
the similarity transformation and whose associated discrimination probability P D is at least as high as Pd- Thus, the 
premise that there is no Ed which gives the highest value of Pd and has the specified symmetry is false. 

It follows from Eq. (3.13), and from the orthogonality of the |7 r ), that the eigenvalues of E D ve are simply Pd/N\c t \ 2 . 
The optimum symmetric operator E° D pt is simply that whose maximum eigenvalue is 1. The desired least upper bound 
on Pd is then given by 

PD<Axmin| Cl .| 2 . (3.15) 

The bound here is clearly less than 1 unless all \c r \ 2 are equal to N , in which case the \ijjj) are orthogonal. 
Although this inequality gives the analytic maximum discrimination probability, in practice it may be necessary to 
employ computational techniques to determine the smallest of the |c r | 2 , as we will see in the next section. To find 
the maximum value of Pd for a specific set of states, it is useful to have an expression for the \c r \ 2 which explicitly 
exhibits their dependence upon the states. Using Eq. (3.8), we find that 

3,0' 

The simplest set of linearly-independent symmetric states comprises just two states. The problem of finding the 
maximum value of Pd for a pair of states has been solved by Ivanovic M , Dieks H and Peres B for equal a priori 
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probabilities and generalised by Jaeger and Shimony |12j] to the case of unequal probabilities. Denoting the two states 
by \i/)±), the Ivanovic-Peres-Dieks limit for the probability of error- free state discrimination is 

iW = i-|^+hMI. (3.i7) 

It is interesting to see how this limit arises as a special case of the bound in (3.15). Up to an irrelevant phase difference, 
the states \tp±) may be represented as 

\ip± >= cos6»|+ > ±sin0|- >, (3.18) 

where the angle 6 lies in the range [0, 7r/4] and the states |± > constitute an orthogonal basis for the space spanned 
by \ip± >■ The system may be represented as a spin-1/2 particle, and |± > taken to be the eigenstates of a z with 
eigenvalues ±1. Note that a z is the unitary operator U relating the states we aim to distinguish, since \ ip± >— ff z \4>^ >■ 
We find that the corresponding reciprocal states are 

\i>± >= sin6>|+ > ±cos0|- > • (3.19) 

The expansion coefficents c± are given by c + = cos# and c_ = sin6*. Within the specified range of 6, |c + | 2 >|c_| 2 , so 
that the maximum value of Pd is 2 1 c | 2 , which is easily seen to be equal to Pidp- 



IV. SYMMETRIC COHERENT STATES 

In this section, we apply the bound (3.15) on Pd to the problem of discriminating between symmetric coherent 
states of the harmonic oscillator or mode of a boson field. These states are 

i„,2 °° q," 

n=0 

where j = 0, N — 1 and aj = ae - " 1 , where a = ao may be any complex number. The \n) are the usual boson 
number states. The magnitudes of the complex arguments aj are all equal to \a\. However their phases are distributed 
around the circle at regular intervals of 2ir/N. Let us denote by P-h the projector onto H, the subspace spanned by 
the \a.j). The unitary transformation U which maps each state onto its successor is 

U = P H e^P H , (4.2) 

where n is the boson number operator. The quantities of interest if we wish to determine the maximum of Pd arc 
the square-moduli of the c r . One can show using Eq. (3.16) that 

| Cr | 2 =i^ e ^ e W 2 (e^-D. (4.3) 

3 

Unfortunately, this summation seems to resist significant simplification, and in general must be carried out numerically. 
A further complication arises if we wish to determine the maximum value of Pd , which entails finding the smallest 
of the \c r \ 2 . This is the fact that for general N, none of the |c r | 2 remains the smallest for all values of \a\ 2 . This can 
be seen in figure 1, which shows the variation of the |c r | 2 as functions of \a\ 2 for N — 10. The behaviour seen here 
for N = 10 is typical of what happens for all N except for N = 2. In this simplest case, \cq | 2 = (l + e- 2 H 2 )/2 and 
| ci | 2 = (1 — e~ 2 l"l )/2, so that |co| 2 >|ci| 2 for all \a\ 2 . For general N, each \c r \ 2 is less than the others for some range 
of \a\ 2 . At \a\ 2 = 0, we find that | Co | 2 = 1 and all of the other c r are zero. Our numerical results for various values 
of N indicate that, as \a\ 2 increases, the smallest of the |c r | 2 is successively |cat— i | 2 , then |ctv— 2 1 2 and so on until it 
is I Co 1 2 then the cycle repeats itself indefinitely. It is evident from the figure that the point at which the minimum 
coefficient changes to a new one occurs when the derivative of the latter is zero. This can be understood when we 
observe that the |c r | 2 obey the relation 

^|v = l^-i| 2 -M 2 . (4.4) 
d(\a\ 2 ) 

It follows that when the derivative of |c r | 2 with respect to \a\ 2 is zero, we have |c r | 2 = |c r _i| 2 . This is the point at 
which these functions cross and thus the smallest function ceases to be |c r | 2 and becomes |c r _i| 2 . As |a| 2 -^oo, the 
|c r | 2 tend to 1/N. Here, the overlaps between the states becomes indefinitely small, and the maximum discrimination 
probability approaches 1. 

The maximum discrimination probability for N — 10 is shown in figure 2. We see that it is an increasing function 
of \a\ 2 , although its derivative is discontinuous whenever a new |c r | 2 becomes the smallest. 
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V. DISCUSSION 



We have obtained the least upper bound on the unambiguous discrimination probability for linearly-independent 
symmetric states with equal a priori probabilities. The corresponding detection operators Epj are the simplest 
possible ones. The only free parameters in these operators for general states are the conditional probabilities Pj. 
For equally-probable symmetric states, the maximum value of Pd is obtained when all Pj are equal, and the one 
remaining free parameter is set by X+(E D pt ) = 1, which is a necessary condition for the maximum value of Pd for 
any linearly-independent set. We thus find that for equally-probable symmetric states, the solution for the maximum 
of Pd is the simplest N state solution. This appears also to be the case for optimisation problems in other areas of 
quantum detection theory, in particular the determination of the unconstrained maximum discrimination probability 
[pj] . p~^| . However, one issue we have not addressed here is whether or not it is only for equiprobable symmetric states 
that setting all Pj to the same value gives the optimum measurement. 

We have also examined the behaviour of the maximum P$ for N symmetric coherent states, in particular its 
dependence on the parameter |a| 2 , which has many interesting features for N > 2. Most notably, the derivative of the 
maximum discrimination probability is not continuous, owing to the fact that no single |c r | 2 is smallest for all values 
of \a\ 2 . While this makes perfect sense from a mathematical point of view, the physical reasons for this pheneomenon 
are by no means obvious and further work may clarify the matter. 
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APPENDIX: CONVEXITY OF THE MAXIMUM EIGENVALUE 

Consider a set of Hermitian operators Ei on a finite-dimensional Hilbert space H., where I = 0, ...,n — 1. Let us 
denote by X+(Ei) the maximum eigenvalue of E\. Then, as we show here, A+ is convex, that is, for any real, positive 
constants a;, 

Z>A+(£i)- (A.l) 
i i 

This admits the following simple proof. Let us write S r = X)I=o a l^l ano - ^ \ s r) be any eigenstate of S r corresponding 
to the maximum eigenvalue X + (S r ). Then, 

A+(SV+i) = (s r+ i\S r + a r E r \s r+ i) 

<X + (S r ) + a r X+{E r ), (A.2) 

where the equality is satisfied only if S r and E r have a simultaneous eigenvector corresponding to the maximum 
eigenvalues of both of these operators. From the definition of S r , we see that 

n— 1 n—1 
a X + (E Q ) + ^[A+(S i+1 ) - X+(Si)} = A+(^ ai E). (A.3) 

1=1 t=0 

Rearranging the inequality in (A.2) gives 

X+(S l+1 ) - X+(§i)<aiX+(Ei). (A.4) 
Substituting this inequality into (A.3) immediately gives (A.l), completing the proof. 
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FIG. 2. Maximum value of the probability Pd of distinguishing between 10 symmetric coherent states as a function of \a[ 
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